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The production of a kaon pair (K+K~ and JsTgisTg) in two-photon processes has been measured using the Belle detector at the 
electron-positron collider KEKB. We have obtained the invariant-mass distribution of the processes between 1.3 and 2.3 GeV. A 
broad bump structure around 1.7 GeV in the K^Kg channel is confirmed, while, in the K + K~ channel, a bump structure is 
seen near 1.9 GeV. The angular distribution of the final-state kaons is analyzed to explore the spin(-helicity) structure in each 
invariant-mass region for the two processes. Contributions from existing and postulated resonances are discussed. 



1 Introduction 

A high-luminosity electron-positron collider is a good 
place to study meson resonances produced by two-photon 
collisions. The well established two-photon resonance in 
kaon-pair final states is the f 2(1525) mesonp'B, which is 
classified as an almost pure ss meson. The L3 Exper- 
iment at LEP has also reported a resonance-like peak 
around 1750 MeV in the K%K% final state I. 

Meson spectroscopy in the 1.5 - 2.0 GeV region is im- 
portant since glueballs are expected to be found in this 
mass range. Glueball candidates around 1.7 and 2.2 GeV 
have been reported in Jj ip — > -yKK decays by the BES 
Experiment at BEPCH. Glueball searches are compli- 
cated by excitations of qq mesons that also populate the 
same region. Since gluons do not couple to photons, the 
two-photon partial decay widths (T 77 ) of pure glueball 
states are expected to be very small. Two-photon pro- 
cesses, therefore, play an important role in the search 
for glueballs. Indeed, null results from the CLEO Ex- 
periment at CESRu lend some support for the glueball 
interpretation of the /j(2220), while the L3 result sug- 
gests the 1700 MeV region should be studied in further 
detail. 



2 Experiment 

The_analysis is based on data taken by the Belle detec- 
tor □ at the KEKB til asymmetric e + e _ collider in the 
period from October, 1999 to December, 2000. The in- 
tegrated luminosity of the data set analyzed for K + K~ 
final startes is 5.86 ib" 1 and 10.6 ftr 1 for K° S K° S final 
states. Since the beam-energy dependence of two-photon 
processes is very small, we combine data taken at the 
Y(4S) resonance (y/s = 10.58 GeV) with data collected 
at values 50 MeV and 60 MeV lower. 

The Belle detector I is a general purpose detector 
which includes a 1.5 T superconducting solenoid magnet 
that surrounds the KEKB crossing point. Charged tracks 



are reconstructed using an 50-layer Central Drift Cham- 
ber (CDC) and a 3-layer double-sided Silicon Vertex De- 
tector (SVD). Particle identification is accomplished by 
combining the responses from a Silica Aerogel Cerenkov 
Counters (ACC) and a Time of Flight Counter system 
(TOF) with specific ionization (dE/dx) measurements in 
the CDC. A Csl Electromagnetic Calorimeter (ECL) lo- 
cated inside the solenoid coil is used for photon detection 
and electron identification. Almost all the signal events 
used in the analysis were triggered by requiring two or 
more tracks in the CDC. 

3 K+K- Final States 

3.1 Event Selection 

Events with only one pair of charged particles produced 
in two-photon processes (77 — > X + X~) are selected 
with the following criteria: the scalar sum of track mo- 
menta (X) m an event is required to be smaller 
than 6 GeV/c, and the sum of the calorimeter ener- 
gies in an event less than 6 GeV. The event is required 
to have only one positively charged track and only one 
negatively charged track, where each satisfies the con- 
ditions: p t > 0.4 GeV/c, \dr\ < 1 cm, \dz\ < 2 cm, 
—0.34 < cos 6 < 0.82, where p t is a transverse mo- 
mentum of a track with respect to the positron beam 
axis, and dr and dz are r and z coordinates, respectively. 
All of the above values are measured in the laboratory 
frame. A cut on p t balance in the e + e~ CM frame, 
|X)Pt| < 0.1 GeV/c, is applied to select exclusive-two- 
track events from quasi-real two-photon collisions. 

We make particle identification cuts for the two 
tracks to select A" + A"~ events. Since the sample included 
a large fraction of 77 — > e + e~ events, we reject them by 
requiring E/p < 0.8 for the two tracks, where E/p is 
the ratio of the energy deposit on ECL to the momen- 
tum. Charged kaons are selected using TOF and ACC 
information with the criteria that the likelihood ratios for 
K/tt and K/p separation are larger than 0.8. In addition, 
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Figure 1: The invariant mass distribution for 77 — ► K + K~ events. 

dE/dx information from CDC is also used for identifying 
kaons in the lowest invariant mass region. 

3.2 Invariant Mass Distribution 

After the application of all the selection criteria, 3674 
events remain. Figure [l] shows the K+K~ invariant 
mass distribution for the selected samples. The peak at 
1.48-1.56 GeV is from .^(1525) -> K+K~. In addition, 
there is a broad bump in the 1.7 — 2.1 GeV region. No 
significant enhancement around 2.23 GeV is observed. 

The mass region near the /j(2220) is shown in detail 
in Figure |[ The signal region is chosen to be within 2220 
MeV and 2260 MeV. To obtain the background shape, 
we fit the K + K~ mass distribution as a second-order 
polynomial from the sideband regions. We observed 48 
events in the signal region and the expected background 
is 49.2. Using a Poisson distribution with background, 
we obtain a 95% C.L. upper limit 

r 77 (/j(2220)) x B(/j(2220) -> K+R-) < 1.5eV 

under the assumption of (J, A) = (2, 2). 

3.3 Analysis of the Angular Distribution 

We made a model-independent partial-wave analysis for 
the angular distribution of the K + K~ events. For the 
total angular momentum (J) of the 77 system, we as- 
sumed that partial waves with J > 4 were negligible. 
The odd- J waves are forbidden by the symmetric nature 
of the initial state and the parity conservation. In the 77 
system, the differential cross section can be written as 

% = \H(J = 0, A = 0) + H{ J = 2, A = 0)| 2 
dil 

+ \H(J = 2,X = 2)\ 2 (1) 
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Figure 2: K + K~ mass distribution (GeV) observed in data near 
the /j(2220) mass. The background shape (dashed line) is assumed 
to be a second-order polynomial. The solid line is the sum of a fit 
to the background and the signal line shape for central values of 
the resonance parameters, Mfj = 2.231 GeV and Yjj =23 MeV, 
corresponding to the observed 95% C.L. upper limit of 14.7 excess 
events. 

where the TVs are the partial-wave amplitudes for the J 
and A (helicity) states. 

Equation (Q) includes four real-number parameters, 
the sizes of the three partial wave amplitudes and one 
relative phase between the two A = waves. Expressing 
Eq. (^) with d-functions, the angular distribution can be 
reduced to 

^^ = A- + 3;|^ (cosr)| 2 + z|^ (cosr)| 2 (2) 

Here the three parameters X , y and Z can be written 
in terms of the three partial-wave amplitudes and one 
relative phase, as explicitly presented in the ReferenceQ. 
However, in Eq. (||), we still cannot determine y and Z 
independently, since the two d-functions have a similar 
shape for |cos(9*| < 0.5. Therefore, we introduced two 
new parameters, Z' = y + Z and / = y/Z, and as- 
sumed / = 0.5. The fitting results for X and Z' are not 
sensitive to the choice of the / value. The parameter X 
contains the J = cross section dominantly, whereas the 
parameter Z' is dominated by the J = 2 contributions, 
if the contribution of the interference effect between the 
two A = states are sufficiently small. 

We have fit the X and Z' values at different W points 
independently. Figure || shows the W dependences of 
the parameters (we give the results in Z'/5 to normalize 
them to the contribution in the total cross section). It 
is clear that the two peak structures from /2(1525) and 
of 1.9 GeV region are dominated by J — 2 components. 
The J = component is consistent with zero except in 
the highest W region. 
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Figure 3: The contribution to the total cross section from the com- 
ponents X (dominated by J = 0, open circles) and 2'/5 (domi- 
nated by J = 2, asterisk marks). 

4 K° S K° Final States 

4-1 Event Selection 

Only events where both Kg's decay to ir + ir~ are consid- 
ered. To select events in which the Kg's were produced 
in two-photon collisions, we require that: the scalar sum 
of track momenta (J^ in an event is smaller than 
6 GeV; the sum of calorimeter energies in an event is 
smaller than 6 GeV; there are exactly four reconstructed 
charged tracks with a net charge of zero; and the total 
momentum imbalance in the transverse plane satisfies 
p t = lEPtl 2 < 0.1 GeV 2 , to select events with small 
incident photon virtuality. 

The charged tracks are assumed to be pions. A pair 
of oppositely charged tracks are combined to form a Kg 
candidate if their reconstructed vertex is displaced from 
the primary interaction vertex. The tracks are then re- 
fit with the constraint that they come from the recon- 
structed vertex. They are kept if the resulting invariant 
mass for the pair is within 10 MeV of M K oE and the com- 
bined momentum vector points back to the interaction 
vertex. After these selection cuts, the mass resolution is 
obtained to be 2.6 MeV. The pairs are again fit, con- 
straining the mass of each Kg candidate to the nominal 
PDG value!. 

Since the two K s 's are produced back-to-back in the 
transverse plane, the angle between the flight directions 
of the two K s candidates in this plane is required to be 
larger than 160°. A legoplot of the two unconstrained 
pion-pair masses is shown in Figure || There is a strong 
enhancement at (m K o , m K o ) and the background frac- 
tion is very low. After applying all the selection cuts, 
1511 events are found in the data sample. 




Figure 4: m w + n — of one Kg candidate versus m n + w - of the other. 

There is a strong enhancement near the (m K o , m K o ) point over a 

S 3 

very small background. 
4-2 Analysis 

The resulting K S K S invariant mass spectrum is shown 
in Figure |[ The spectrum obtained-is similar to the 
one reported by the L3 Experiment □. The spectrum 
is dominated by the /2(1525) resonance and a clear 
enhancement is visible in the 1750 MeV region. The 
f% ( 1270) — a\( 1320) region shows the destructive interfer- 
ence expected in the K S K S final state □. The event rate 
in this region is suppressed by the detection efficiency 

The K S K S invariant mass spectrum is fitted by min- 
imizing a x 2 function for the expected value in the ith 
bin. 

Ek = & • (Gi + BW} 1] + BW} 2) + Bi) . (3) 

where £j is the detection efficiency parameterized from 
the Monte Carlo. The small bump in the / 2 (1270) - 
a2(1320) region is described by a Gaussian function d 
with its mean and sigma as free parameters. The peaks 
at /2(1525) and 1750 MeV regions are described by Breit- 
Wigner functions BWi with their masses and widths 
taken as free parameters. The background is assumed 
to be a constant (Bi) in the fit. The result of the fit 
is illustrated by the curve in Figure |5[ The parameters 
obtained are summarized in Table 

4.3 The f j (2220) Mass Region 

No enhancement at the /,/(2220) mass is observed. Wc 
assume a mass of 2230 MeV and a total width of 20 MeV. 
The signal region is chosen to be ±40 MeV around 2230 
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Figure 5: The K ( ^K^ invariant mass spectrum for the 77 
process. The solid line is a fit with two Breit-Wigner functions 
for the / '2(1525) and 1750 MeV region, and a Gaussian for the 
/2(1270) — a2(1320) region plus a constant background. See the 
text for details. 



Table 1: Results of the fit to the K°K° mass spectrum. 





#(1525) 


1750 MeV 


Mass (MeV) 


1532 ± 4 


1768 ±9.6 


Width (MeV) 


64 ±6.8 


323 ± 29 


No. of Events 


414 ±36 


967 ± 72 
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Figure 6: The contribution to the total cross section from the com- 
ponents (a) Z' /5, dominated by spin-2 component, and (b) X, 
dominated by spin-0 component. 



large spin-0 component in the 1750 MeV region while 
some spin-2 structure appears around 1.8 — 2.0 GeV. 

5 Summary 

We have studied the reactions 77 — * K + K~ and 77 — ► 
KgKg using the large data samples collected by the Belle 
experiment at KEKB. A prominent /2(1525) resonance is 
observed in both channels. A broad structure in the 1.7— 
2.1 GeV region is found in K + K~ and an enhancement 
around 1750 MeV is observed in KgKg. A partial wave 
analysis has been performed to explore these resonance 
structures. Upper limits for /j(2220) are also obtained 
in both channels, respectively. 
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sis, described in Sec. |3.3| , to the KgK§ 
angular distribution of each W bin is fit by Eq. 
ure|^ shows the W dependences of the fitted results, Z'/5 
and X , which correspond to the contribution of J = 2 
and J = components in the total cross section, respec- 
tively. It is clear that the peak of /2(1525) is dominated 
by spin-2 component as expected. However, there is a 
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